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Abstract 

We present a stochastic model of dynamically interacting grains in one dimen- 
sion, in the presence of a low vibrational intensity, to investigate the effect of shape 
on the statics and dynamics of the compaction process. Regularity and irregularity 
in grain shapes are shown to be centrally important in determining the statics of 
close-packing states, as well as the nature of zero- and low-temperature dynamics 
in this columnar model. 
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That shape matters crucially in the efficient packing of objects is everyday experience. 
Coffee grains can be shaken into a closely packed state far more readily than marbles, 
where necks and pores remain even in the most ordered state. This Letter is devoted to 
a simple stochastic dynamical model in one dimension, demonstrating the effect of shape 
in closely packed systems, such as glasses or densely packed granular media (see (T] for 
a review, and for recent investigations of lattice models of granular systems). A 
more general version of the model, where lattice sites can be empty, has been described 
elsewhere Oil]. The present model, and some of its generalisations, will be investigated 
in more detail in 

The model 

Our model consists in a column of sites, each of which is occupied by a grain. Grains 
are indexed by their depth n = 1, . . . ,N, measured from the top of the column. Grains are 
assumed to be anisotropic in shape. They take, for simplicity, two possible orientations, 
referred to as ordered and disordered. We define orientation variables by setting (T„ = +1 
if grain number n is ordered, and (T„ = — 1 if grain number n is disordered. A configuration 
of the system is uniquely defined by the A^ orientation variables {cr„}. While ordered 
grains are perfectly packed, disordered grains are imperfectly packed. Each disordered 
grain leaves a void space e on the site it inhabits. 

Non-trivial dynamical interactions between the grains are such as to minimise the void 
space locally. These lead to frustration between competing grain orientations which seek 
to minimise interstitial voids. Under the influence of a dimensionless vibration intensity F, 
this system undergoes cooperative reorganisations reminiscent of glassy dynamics [U]. 
From a phenomenological viewpoint, we model the above by a stochastic dynamics, 
defined by the following transition rates per unit time: 

( w„(cr„ = + cr„ = -) = exp {-n/^dyn - hn/T) , 

(1) 

[ Wn(0"n = > Cr„ = +) = exp (-n/^dyn + hn/T) . 

The dynamical rules (Q) fully define the model. In these formulas, F is a dimensionless 
measure of the vibration intensity, which will also be referred to as temperature. The 
dynamical length ^dyn is a phenomenological parameter which controls the spatial depen- 
dence of dynamical behaviour (In earlier work |^ ^dyn was shown to determine the 
extent to which order propagates down the column in the glassy regime, a perspective 
that is maintained here). Finally, and most importantly, the ordering field hn acting on 
the orientation cr^ of grain number n is compacting (see later): it reads 

hn = em~ -m^, (2) 

where m+ (resp. m~) is the number of ordered (resp. disordered) grains above grain 
number n: 

n—l n— 1 

= 15 E(l + "ni.- = -Y,{l-ak), (3) 

^ k=i ^ k=i 
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so that m+ + m~ = n. Hence represents the excess void space j7] of the system. 
Equation (j2)) shows that a transition from an ordered to a disordered state for grain 
number n is hindered hj the number of voids that are already above it, i.e., represents a 
compacting dynamics. Our model is thus fully directed in time and space: the orientation 
of a given grain only influences the grains below it and at later times. The free surface 
of such a system under vibration could be expected to have maximal mobility; the free 
volume generated by vibration fSj would decrease with increasing depth in the column. 

In the limit of a vanishing vibration intensity (F 0), the probabilistic rules (0) 
become the following deterministic updating formula: 

cr„, sign hn = sign(£m^ - "^n)) (4) 
provided hn 7^ (see below). 

Statics 

Ground states are defined as being the configurations such that the static counterpart 
of dH): 

cr„ = sign hn = sign(£:m^ - m+), (5) 

holds for all grains n = 1, . . . , N. In other words, in a ground state each orientation cr„, 
is aligned along its local field hn- 

We mention for completeness that the case e < is a generalisation of earlier work 0, 
with a complete absence of frustration and a single ground state of ordered grains. In 
the present situation {e > 0), a. rich ground-state structure is achieved, because of frus- 
tration jUj, whose nature depends on whether e is rational or irrational. The rotation 
number 

Q = e/{e + l) (6) 

fixes the proportions of ordered and disordered grains in the ground states: /+ = Q, 
/_ = 1 — f2, and their difference = f+ — f- = 2Q — 1 can be seen as a spatially 
averaged 'magnetisation'. 

• For irrational e, (j2I) implies that all the local fields hn are non-zero. A unique quasiperi- 
odic ground state is thus generated 0. It can be constructed by the cut-and-project 
method, along the lines of a geometrical approach developed for quasicrystals JU]- The 
local fields /i„ lie in a bounded interval —l<hn<£- 

• For rational e = p/q, with p and q mutual primes, Q = p/{p + q), and some of the 
hn can vanish. The equation hn = means that grain number n has a perfectly packed 
column above it, so that it is free to choose its orientation. For e = 1/2, for example, 
one can visualise that each disordered grain carries a void half its size, so that units 

of perfect packing must be permutations of the triad H , where the two half voids 

from each of the — grains are filled by the + grain. The dynamics selects two of these 

patterns, H and — I — . More generally, orientational indeterminacy occurs at points 

of perfect packing such that n is a multiple of the period p + q. Each ground state is 
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a random sequence of two patterns of length p + q, each containing p ordered and q 
disordered grains. The model therefore has a zero-temperature configurational entropy 
(or ground-state entropy) per grain: 



E = \n2/{p + q). 



(7) 



Associating (ir)regular grains with (ir)rational void spaces, this has the appealing 
physical interpretation that irregularities in grain shapes lead to a unique state of close 
packing (such that all jagged edges are well meshed together), while regular grains have 
a huge degeneracy of such states (as in the fabled greengrocer's problem [H]). 

Zero-temperature dynamics 

We now turn to the zero-temperature dynamics (@)), starting with a disordered initial 
configuration. 

• For irrational e, the application of this dynamics causes the quasiperiodic ground 
state to be recovered downwards from the top of the column. The depth of the ordered 
boundary layer grows ballistically with time: 



The velocity V{e) = V{l/6) varies smoothly with e, and diverges as V{e) ~ e for 
£ ^ 1 ini- The rest of the system remains in its disordered initial state. When L{t) 
becomes comparable with ,^dyn, the effects of the free surface begin to be damped. In 
particular for t ^ ^dyn/V{e) we recover the logarithmic coarsening law L(t) ^ ^dynlnt, 
observed in related work ^] to model the slow dynamical relaxation of vibrated 
sand [THj . 

• For rational e, as mentioned above, the local field may vanish. The corresponding 
orientation is updated according to (T„ — > ±1 with probability 1/2, leading to a dynamics 
which is stochastic even at zero temperature. Here, even the behaviour well within the 
boundary layer ,^dyn contains many intriguing features, while the dynamics for n ^ ^dyn 
again takes place on a logarithmic scale We therefore focus on the limit ^dyn = oo. 
The main result is that zero-temperature dynamics does not drive the system to any of 
its degenerate ground states. The system instead shows a fast relaxation to a non-trivial 
steady state, independent of its initial condition. The local fields hn have unbounded 
fluctuations in this steady state, which are reminiscent of the density fluctuations about 
the mean packing fraction observed in granular systems [13 Ej , especially around the 
so-called random close packing density [7| ? which is the highest density achievable in 
practice by extensive dynamical processes. 

Figure ^ demonstrates the anomalous roughening law 



The local fields hn are approximately Gaussianly distributed. The following scaling 
argument explains the observed roughening exponent 2/3. Let /i„ be the position of a 





= {hi) ^ An'/-' 



{A ^ 0.83). 



(9) 
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Figure 1: Log-log plot of = (/i^) against depth n, for zero-temperature dynamics with 
e = 1. Full line: numerical data. Dashed line: fit of asymptotic behaviour, leading to Q- 

fictitious random walker at time n. The noise in this random walk originates in the sites 
m < n where the local field vanishes. It is therefore proportional to Z]m=\ Prob{/im = 
0}, hence the consistency condition ~ I]m=i l/W^m, yielding the power law Q. A 
more detailed derivation will be given in 0. Roughly speaking, the walker obeys a 
diffusion equation, with a diffusion constant scaling as the inverse distance to its starting 
point. This is reminiscent of the domain-growth mechanism in the low-temperature 
coarsening regime of the Ising chain with Kawasaki dynamics the power law L{t) ~ 
t^/^ for the mean domain size (analogous to Wn ~ n^/^) can be understood from the 
picture of diffusing domains, whose diffusion constant scales as the inverse of their length. 
The anomalous roughening law Q is the most central feature of the zero-temperature 
steady state observed for rational e. 

If the grain orientations were statistically independent, i.e., uncorrelated, one would 
have the simple result (/i^) = (e^ -|- l)n/2, while Q implies that {hi) grows much more 
slowly than n. Orientational displacements are thus fully anticorrelated. More precisely, 
Figure 121 shows that the orientation correlations Cm,n = {(^m<^n) scale as 



The scaling function F is even, positive, and it obeys F{x) dx = 1, expressing that 
spin fluctuations are asymptotically totally screened (i.e., fully anticorrelated). This is 
similar to the bridge collapse seen in displacement-displacement correlations of strongly 
compacting grains p^; grain orientational displacements in the direction of vibration 





(10) 
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Figure 2: Scaling plot of the orientation correlation function Cm,n for n 7^ m in the zero- 
temperature steady state with e = 1, demonstrating the validity of and showing a plot of 
(minus) the scaling function F. 



were there seen to be strongly anticorrelated in jammed regions, as each grain tried to 
collapse into the void space trapped by its neighbours. We remark that temporal anti- 
correlations have also been observed in recent experiments investigating the properties 
of cages near the colloidal glass transition Interestingly, correlations transverse to 
the shaking direction were found to be rather small ^3], thus, in self-consistency terms 
justifying the choice of a column model in the present case. 

Figure El shows a plot of the normalised probabilities 2^p(C) of all the configurations C 
in the zero-temperature steady state of a system with = 12 grains and e = 1, directly 
measured in a very long simulation. The probabilities are plotted against the 2^^ = 4096 
configurations, sorted lexicographically (read down the column). This plot exhibits a 
rugged structure on this microscopic scale: some configurations are clearly visited far 
more often than others. We suggest that this behaviour is generic: i.e., the dynamics of 
compaction in the jammed state leads to a microscopic sampling of configuration space 
which is highly non-uniform. In spite of this fine structure, the steady-state entropy 

^ = -^p(C)lnp(C), (11) 
c 

is not far from the estimate 5'fiat = In 2 along the lines of Edwards' flatness hypothe- 
sis jTHJ. Indeed, for = 12 (data of Figure El), we have 5* = 7.839 against ^flat = 8.318. 
In other words, the entropy reduction jTH] AS = Sfiat — S = 0.479 is small, compared to 
the absolute value of the entropy per grain. 
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Figure 3: Plot of the normalised probabilities 2^ p{C) of all the configurations in the zero- 
temperature steady state for = 12 and e = 1, against the configurations C sorted lexico- 
graphically. 

Low-temperature dynamics 

Low-temperature dynamics induces drastic changes in the case of e irrational. (Since 
zero-temperature dynamics remains stochastic in the rational case, we do not expect 
low-temperature dynamics to introduce qualitative differences there.) For a low but non- 
zero r, there will be a few mistakes, i.e., orientations which are not aligned with their 
local field according to (0). The a priori probability of observing a mistake at site n 
scales as 

nH ~exp(-2|/i„|/r). (12) 

Hence the sites n so that \hn\ ~ F ^ 1 will be preferred nucleation sites for mistakes, and 
thus dominate the low-temperature dynamics. It turns out that those sites are such that 
nQ is close to an integer. The magnitude of the excess void space is least here, leading 
to the least cost for a misalignment. 

Choosing the golden mean e = $ = (-\/5 + l)/2 for concreteness, the preferred nu- 
cleation sites are given by the Fibonacci numbers: n = ^ ^'^/a/S. Then n(Ffc) ~ 
exp(— 2$^/(-\/5 FFfc)). Let us denote the instantaneous position of the uppermost mis- 
take by A/'(t). This depth divides an upper boundary layer, ordered in the quasiperiodic 
ground state, from a disordered lower zone, characterised by the anomalous roughening 
law 0- The mistake is itself advected ballistically with velocity V{^) ~ 2.57, just as 
with zero-temperature dynamics, until another mistake is nucleated above it. 

Figure IH shows a typical sawtooth plot of the instantaneous depth A/'(t) , for a tem- 
perature F = 0.003. The ordering length (TV) is expected to diverge at low temperature. 
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Figure 4: Plot of the instantaneous depth M{t) of the ordered layer, for e = $ (the golden 
mean) and T = 0.003. Dashed lines: leading nucleation sites given by Fibonacci numbers 
(bottom to top: Fn = 89, Fu = 144, F13 = 233). 

as mistakes become more and more rare. From a more quantitative viewpoint, the most 
active Fibonacci site F^, is such that the nucleation time l/Il{Fk) is comparable to the 
advection time to the next nucleation site -Ffc+i- This yields 

(Ar)~i/(r|inr|). (13) 

The ordering length therefore diverges linearly at low temperature, up to a logarithmic 
correction. A similar law is predicted for all the irrational values of e with typical Dio- 
phantine properties [5j. This mechanism, where mistakes appear at preferred nucleation 
sites, is reminiscent of the phenomenon of hierarchical melting observed in incommensu- 
rate modulated solids PU] . 

Both (A/") and .^dyn retain the flavour of a boundary layer separating order from 
disorder. Within each of these boundary layers, the relaxation is fast, and based on 
single-particle relaxation, i.e., individual particles attaining their positions of optimal 
local packing ^21 ^M- The slow dynamics of cooperative relaxation only sets in for 
lengths beyond these, when the lengths over which packing needs to be optimised become 
non-local. This in turn leads, as in reality [7], to hysteresis, i.e., a dependence on the 
initial state of the packing. 

Remarkably, all of these features were obtained at a qualitative level in the glassy 
regime of a much simpler model [^j. On the one hand, this allows us to speculate 
that the shape-dependent aging phenomena seen there could be retrieved here, i.e., that 
conventional aging phenomena would only be seen for irregular grains (irrational e). 
On the other hand, it is tempting to ask if the directional causality of the dynamical 
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interactions present in this model and the earher one could be responsible for their 
qualitative similarity, and thus be a necessary ingredient for modelling glassiness. 
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